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ABSTRACT 

We investigate the effect of spiral structure on the Galactic disk as viewed by pencil beams centered 
on the Sun, relevant to upcoming surveys such as ARGOS, SEGUE, and GAIA. We create synthetic 
Galactic maps which we call Pencil Beam Maps (PBMs) of the following observables: line-of-sight 
velocities, the corresponding velocity dispersion, and the stellar number density that are functions of 
distance from the observer. We show that such maps can be used to infer spiral structure parameters, 
such as pattern speed, solar phase angle, and number of arms. The mean line-of-sight velocity and 
velocity dispersion are affected by up to ~ 35 km/s which is well within the detectable limit for 
forthcoming radial velocity surveys. One can measure the pattern speed by searching for imprints 
of resonances. In the case of a two-armed spiral structure it can be inferred from the radius of a 
high velocity dispersion ring situated at the 2:1 ILR. This information, however, must be combined 
with information related to the velocities and stellar number density in order to distinguish from a 
four-armed structure. If the pattern speed is such that the 2:1 ILR is hidden inside the Galactic bulge 
the 2:1 OLR will be present in the outer Galaxy and thus can equivalently be used to estimate the 
pattern speed. Once the pattern speed is known the solar angle can be estimated from the line-of-sight 
velocities and the number density PBMs. Forthcoming radial velocity surveys are likely to provide 
powerful constraints of the structure of the Milky Way disk. 
Subject headings: stellar dynamics, spiral structure 



1. INTRODUCTION 

It has been well established by now that the Milky 
way is not axisymmetric with both a central bar and 
spiral structure perturbing its disk. Due to our location 
in the Galactic plane both spiral and bar structure is 
impossible to observe directly. Galactic bar parameters 
such as orientation and pattern speed have been inferred 
indirectly f rom both asymmetries around the Ga lactic 
center fe.g.. lBlitz and SpergellH99l iWeinber J[l99l and 
its effect on the local vel ocity distribu t ion o f old stars , 
i.e., the Hercules str eam (|Dehnenlll99l l2000t IPu"xll2001l : 
iMinchev et~alll2007h . 

Spiral structure parameters, however, are much 
more unce r tain. Current spiral density wave models 
(iFuxl 120011: ILepine et alJ 120011 : iDe Simone et alJ 12004 
iQuillen and Minchevl l2005f ) strongly disagree on the 
strength of the spiral structure, the number of arms, and 
the pattern speed. These models differ in their predic- 
tions of the induced velocity streaming at different angu- 
lar positions in the Galaxy. For example, a four-armed 
density wave with velocity perturbations of ~ 20 km/s 
will exhibit rapidly varying radial and tangential velocity 
components with azimuth across distances of a few kpc, 
and we could expect to detect ^20 — 50 km/s variations 
in the mean line-of-sight stellar velocity as a function 
of the distance from the Sun. However, the strength of 
the spiral arm perturbation re mains controvers i al. Ba sed 
on the OGLE number counts, iPaczvnski et al.l (|1994l ) es- 
timated that the Sagittarius-Carina arm has a factor 
of two increase in density compared to the underlying 
disk. This model is inconsistent with COBE studies 
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which find a much smaller contrast (~ 15%) and show 
that the Perseus and Scu-C ru arms are more dominant 
(jDrimmel and Sp"e rgel 200j]). 

HI, CO, Cepheid, and far-infrared observations sug- 
gest that the Galactic disk contains a four-armed 
tightly wound stru c ture. On the other hand, 
iDrimmel and Sperg el (2001) have shown that the near- 
infrared observatio ns are consist e nt wit h a dominant two- 
armed structure. ILepine et al.l (|2001h suggest that lo- 
cally the Milky Way can be modeled by the superpo- 
sition of a two- and four-armed structure moving at the 
same pattern speed. By studying the nearby spiral arms, 
iNaoz and~S haviv (2007) find that the Sagittarius-Carina 
arm is a superposition of two features, moving at differ- 
ent pattern speeds. The effect of a two- and four-armed 
structure, moving at different angular velocities, on the 
velocity disper sion of a galactic disk has b een explored 
numerically by IMinchev and Quillenl (|2006( ) . 

Estimates for the pattern speed of the Milky Way spi- 
ral structure, or equivalently, the Sun's position with re- 
spect to resonances associated with spiral structu re, span 
a larg e range of values. Reviewing previous work. IShavivl 
(2003) finds a clustering of estimates for the pattern 

speed of local spiral structure near O s ~ 20kms~ 1 kpc , 
though ot her studies suggest f l s ~ 13kms _1 kpc . The 
model by ILepine et al.l (|2001l ) places the Sun near the 

corotation resonance Q s ~ 28kms~ 1 kpc 1 ), and was fit 
to Cepheid kinematics. The recent gas dyna mical stud- 
ies (jMartos et al.l 12004 iBissantz et aT1l2003l ) match the 
properties of the gas in nearby arms with a spiral pat- 
tern speed of ~ 20kms _1 kpc _1 . iMartos et alJ (|2004f ) pro- 
pose that a two-armed stellar structure consistent with 
the stellar distribution inferred from COBE could cause 
four-arms in the gas distribution near the Sun. The 
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gas dynamical model proposed by iBissantz et alj (|2003f ) 
with a similar spiral pattern speed matches HI and CO 
kinematics. The pattern speed of a spiral density wave 
can be tightly constrai ned from the location of its res- 
onances. For example, iQuillen and Minchevl (|2005l ) as- 
sociated stellar streams in the solar neighborhood with 
the 4:1 ILR resonance of a two-armed pattern and were 
then able to tightly constrain the pattern speed of the 
driving spiral density wave to within 5%. Independent 
constraints on the pattern speed come from recent sur- 
veys of nearby open clusters (e.g. iDias and LepinefeOOSl ) 
where the older clusters are found to have drifted further 
from their original density wave location. These authors 
concluded that the Sun is located n ear the CR. A solar 
circ le near the CR is also fav ored bv lLepine et alj (|2001f ) 
and lNaoz and Shavivl (|2007f ). 

In this paper we investigate how spiral structure pa- 
rameters can be inferred from velocity and density maps 
resulting from pencil-beam and large-scale surveys of the 
Galaxy. At present the influence of spiral arms on the ob- 
served kinematic properties of the Galactic disk is very 
poorly understood. With the advent of future Galac- 
tic all-sky (GAIA, SEGUE) and pencil-beam (ARGOS, 
BRAVA) radial velocity surveys, large amounts of kine- 
matic data will be collected. The types of dynamical 
constraints made possible with these new data sets is 
not currently known. We address that issue here with 
synthetic models for the purpose of exploring how spiral 
structure might be constraint from these data. 

2. THE SIMULATIONS 

We perform 2D test-particle simulations of an initially 
axisymmetric exponential galactic disk. In order to re- 
produce the observed kinematics of the Galactic disk, 
we use disk parameters consistent with observations (ta- 
ble [it . The reader is referred to iMinchev and Quillenl 
(2007) for a more detailed description of our simulation 
set up. In all of our simulations we start with an ini- 
tially warm disk, i.e., the radial velocity dispersion at ro 
is a u — 0.20«o where vo is the velocity of the local stan- 
dard of rest. The background axisymmetric potential 
due to the disk and halo has the form $o(r) = Vq log(r), 
corresponding to a flat rotation curve. 

We treat the spiral pattern as a small perturbation to 
the axisymmetric model of the galaxy by viewing it as 
a quasi-steady d ensity wave in a ccordance with the Lin- 
Shu hypothesis (jLin et al.l fl969). The spiral wave grav- 
itational potential perturbation is expanded in Fourier 
components as 

$i(r,0,t) = ^e m exp [i(abxr - m(0 - . (1) 

m 

The parameter a is related to the pitch angle of the spiral 
wave, p, as a — mcot(p), negative for trailing spirals 
with rotation counterclockwise, and (r, <j)) are plane polar 
coordinates. The pattern speed is given by fi s and the 
spiral strength by e m . For a two- armed structure the 
m = 2 term dominates. Upon taking the real part of 
equation [T] the perturbation due to the two- armed spiral 
density wave becomes 

^(r,^,*) = e s cos(alnr-2(^-n s t)). (2) 

Integrations are performed forward in time. The per- 
turbation is grown from zero to its maximum strength in 



TABLE 1 
Simulation parameters used 



Parameter 


Symbol 


Value 


Solar neighborhood radius 


ro 


1 


Circular velocity at ro 


v 


1 


Radial velocity dispersion 


c u (ro) 


0.20v 


tr u scale length 


T a 


0.9r 


Disk scale length 


r P 


0.37r 


Spiral strength 


e s 


-0.015 


Pitch angle 


P 


18° 



four rotation periods at tq. In order to improve statis- 
tics, positions and velocities are time averaged for 10 
spiral periods. We distribute particles (stars) between in 
inner and outer galactic radii {ri n ,r ou t) = (0.3ro, 2.0rrj)- 
New particles are added until the final number of outputs 
is 2.5 x 10 6 . In addition, the two-fold symmetry of our 
model galaxy is used to double this number. 

We present our results by changing the spiral pattern 
speed, Q s , and keeping the solar radius fixed at r = 1. 
For a two-armed spiral pattern the primary resonances 
are the 2:1 inner and outer Lindblad resonances (ILR 
and OLR). Those are achieved when Q s /f2rj = 1|k/2« 
0.3, 1.7, respectively, where n is the epicyclic frequency. 
Similarly, the second order resonances are the 4:1 Lind- 
blad resonances (LRs) at f2 s /f2o = 1 + K /2 ~ 0.65, 1.35. 
We examine a region of parameter space for a range of 
pattern speeds placing the SN between the 4:1 LRs. 

3. VARIATION OF GALAXY MORPHOLOGY WITH 
PATTERN SPEED 

Interpretation of line-of-sight velocities with Galactic 
longitude and distance from the Sun is not straightfor- 
ward. To help out we first discuss morphology as seen 
by an outside viewer. 

In figure [1] we present stellar number density contour 
plots for simulations of galactic disks with different pat- 
tern speeds. The background axisymmetric disk is sub- 
tracted to emphasize the spiral structure. The quantity 
plotted is (S — T, ax i)/T, ax i, where £ and ~S a xi are the per- 
turbed and axisymmetric stellar number densities. Con- 
centric circles represent the 2:1 LRs (dashed), the solar 
radius (solid), and the CR (dash-dotted). Darker col- 
ors correspond to lower density. The inner 0.3ro disk is 
not plotted since we do not model the Galactic center. 
Each panel represents a simulation with a distinct pat- 
tern speed, f2 s , and all other parameters kept the same 
(see table [T]). Pattern speeds considered range approx- 
imately between the 4:1 LRs, f2 s = [0.6, 1.3] f2o m units 
of 0.1. The minima of the two-armed spiral potential are 
graphed in each panel as solid curves. Note the crowding 
of resonances as the pattern speed increases. 

In general, changing the solar radius in a simulation 
with the same pattern speed is equivalent to changing the 
pattern speed and keeping the solar radius fixed. How- 
ever, this is exactly true only if the stellar density and ve- 
locity dispersion varied linearly with radius. This is not 
the case in real galaxies; b oth of these are found to v ary 
exponentially with radius. iLewis and Freemanl (|1989h es- 
timated the number density and radial velocity disper- 
sion scale lengths in the Milky Way to be r p = 0.37ro, 
and r a = 0.9ro, respectively. Thus we need to per- 
form different simulation runs when changing the pattern 
speed. 




Fig. 1. — Each panel shows stellar number density contour plots of a simulation with a particular spiral pattern speed, Q s , in units of the 
local circular velocity, Oo- The axisymmetric disk is subtracted to emphasize the perturbation. The quantity plotted is (S — Y, ax i)/H ax i , 
where S and T, ax i are the perturbed and axisymmetric stellar number densities. Contours are spaced linearly with white indicating highest 
values. Dashed circles show 2:1 LRs, the dash-dotted line indicates the CR, and the solid circle represents the solar neighborhood at 
r/ro = 1. The inner 0.3ro is not plotted since we do not simulate the central bar. We consider pattern speeds between the 4:1 ILR and 
the 4:1 OLR. 



Note the disruption of the spirals near the 2:1 LRs 
(dashed circles). A rapid decrease of spiral strength at 
the 2: 1 OLR was also observed by iVorobvov and Theia 
(pool in a galactic disk model consisting of solving nu- 
merically the Bolt zmann moment equ ations. It has also 
been suggested bv lContopoulod (|1985[ ) that strong (non- 
linear) spiral structure cannot extend beyond the 4:1 ILR 
since at that location the stellar orbits are not in phase 
with t he imposed spiral. As pointed out by ISellwoodl 
(1993) , however , this limited extent of the spirals found 
bv lContopoulosI (|1985[ ) is probably related to the restric- 
tive assumptions they make in order to constructing self- 
consistent spiral structure. Another remarkable feature 
in the plots of figure Q] is the overdensity of stellar or- 
bits just outside the 2:1 LRs (where they exist) and near 
the CR (dash-dotted circles). In the case of the CR, the 
enhancement is due to the stable Lagrange points Ln 
associated with this resonance. 

Assuming the primary spiral structure in the solar 
neighborhood is two-armed, the question arises: How can 
we determine any spiral structure parameters, given our 
inconvenient position in the Galaxy? One way to do this 
is by collecting a large number of stars with known veloc- 
ities, distances, etc, and constructing velocity and den- 
sity maps by plotting these versus Galactic longitude, I 
and heliocentric distance, d. In the next section we show 
what such maps would look like when a two-armed spiral 
structure perturbs a stellar disk. We examine different 
pattern speeds and solar positions with respect to a spi- 
ral arm. We call these maps "Pencil Beam Maps" or 
PBMs. 



4. PENCIL BEAM MAPS (PBMS) OF A GALACTIC DISK 
PERTURBED BY A TWO-ARMED SPIRAL STRUCTURE 

To investigate the global structure of the Galaxy we 
require accurate stellar velocities and distances. In a 
pencil-beam spectroscopic survey line-of-sight velocities 
can be measured to great distances. On the other hand, 
proper motions are hard to measure for stars farther than 
about two kpc from the Sun. Thus those cannot be used 
in our investigation. 

For a complete kinematic study accurate distance es- 
timates are also needed. Due to the large distances in- 
volved in a such survey trigonometric parallax measure- 
ments are not possible. Instead, photometric distances 
can be estimated given accurate photometry. This way 
of computing distances, however, is hampered by the 
dust obscuration in the Galactic plane aside from sev- 
eral known windows, e.g., Baade's Window at (I, b) — 
(0.9°,— 4°). Another distance estimator is the use of 
standard candles such as Cepheids, etc. Here we do not 
attempt to model the reddening resulting from dust ex- 
tinction but present an idealized model as a first attempt 
to tackle this problem. A future paper will be dedicated 
to a more detailed modeling. Due to this shortcoming 
our model can be directly applied only to the known low 
extinction Galactic plane windows. Like Baade s win- 
dow, the Scutum window at I — 27° has low extinction 
and we can observe stars at 10 kpc distances towards 
the inner disk. Clump giants of the intermediate-age 
and older population of the disk and thick disk will be 
abundant in these fields. This line-of-sight at I = 27° is 
tangent to the Scu-Cru spiral arm, with an AV extinc- 
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Fig. 2. — Pencil-beam maps of an axisymmetric disk (first row), same with an added two-armed spiral structure perturbation (second 
row), and the sole effect of the spiral structure (third row). The pattern speed is Cl a = 0.7f2o and the observer is at a solar phase angle 
of <po = 20°. Columns from left to right show contour plots of the line-of-sight velocity v^, its standard deviation u Vd , and the number 
density S. All of these are plotted versus the Galactic longitude I (x-axis), and the heliocentric distance d/ro (y-axis). The white curve 
and vertical line are the projection of the solar circle and the Galactic longitude 1 = 0°, respectively. The number density plot in the x-y 
plane of this simulation was shown in figure HI (panel with Q s = 0.7). Clearly, the nonaxisymmetric structure caused by the effect of the 
spiral arms is much better pronounced in the axisymmetric background subtracted values (third row). 



tion of about 3 mag at the distance of the spiral arm 
tangent point (~ 6 kpc). In the Scutum window, the HI 
and Ha profiles clearly show the presence of spiral arms 
(jMadsen and Revnoldsll2005fh 

In figure[2]we present PBMs of the line-of-sight velocity 
Vd (left column), the corresponding velocity dispersion 
a Vd (middle column), and the stellar number density £ 
(right column). This is a simulation of a galactic disk 
perturbed by a two-armed spiral density wave moving 
with fl s = 0.7il . To create the contour plots in this fig- 
ure we bin the disk in Galactic longitude I (x-axis), and 
heliocentric distance d/ro (y-axis) as seen from an ob- 
server at a solar orientation with respect to the concave 
spiral arm of 0o = 20° . This is in contrast to figure Q] 
where we present number density plots of a face-on view. 
The contours in the first row in figure [2] show the results 
of an axisymmetric disk, which is indicated by the sub- 
script "axi" . In the second row the disk is perturbed by 
an m = 2 spiral density wave. The third row in figure 
[2] plots contours of the difference between the perturbed 
and axisymmetric disks for the mean velocity and its 
dispersion: Av d = v d - v d ,axi and Aa Vd = a Vd - a Vd ^ axi . 
The number density, on the other hand, is obtained as 
in figure [TJ AE = (£ — T, axi )/Yj axi We showed a num- 
ber density plot in the x-y plane of this particular sim- 
ulation in figure [T] (panel with fl s — 0.7). In figure [2l 



however, we plot observables from a point of view cen- 
tered on the Sun, as pencil-beam surveys would see the 
Galaxy. The shaded contours are equally spaced with 
darker color corresponding to lower density. On top of 
each panel we show the minimum and maximum contour 
values (in units of vq in the case of the velocities and the 
velocity dispersion). As it is commonly accepted, the 
Galactic longitude is zero in the direction of the Galactic 
center, with the Galactic anticenter at I = ±180°. The 
inner 0.3ro disk has been removed everywhere as in the 
density plots in figure [1] The white curve in each panel 
represents the projection of the solar circle in a PBM and 
the vertical line shows the Galactic longitude I = 0°. 

Note that in the background subtracted PBMs (third 
row) the spiral structure is much better pronounced com- 
pared to the raw values (first row). Knowing the global 
Galactic potential is imperative to extracting informa- 
tion from a PBM. This presents a problem since the 
Milky Way potential is not very well known. We could, 
however, use our model axisymmetric data to subtract 
from the observational data. 

How well do we need to know the axisymmetric struc- 
ture so that spiral features are not wiped out? From the 
third row of figure [2] we can estimate that distance un- 
certainties of ~ 30% do not prevent detection of spiral 
structure. We also need line-of-sight velocity precision 
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Fig. 3. — Pencil-beam maps showing variation of structure in the observables with the change in pattern speed. The solar angle is kept 
fixed at <j>o = 20° everywhere. In contrast to figure [2] here all panels have the axisymmetric background subtracted to reveal the symmetry 
of the spiral residuals. Note for the case of A(T„ d the ring-like shape around the Galactic center for all pattern speeds (middle column). 
These "hot" rings are associated with the 2:1 ILR induced by the spiral density wave. As the 2:1 OLR enters the disk (around Q s = 0.9f2o, 
see figure [TJ strong features in Ao\, d PBMs appear in the outer disk. 



of ~ 20 km/s and the axisymmetric potential must be 
known to within ~ 10%. 

What information about the spiral structure can we in- 
fer from PBMs such as figure [2]? To answer this question 
we need to vary the parameters and look at how struc- 
ture in these PBMs changes. We do this in the following 
sections. 



4.1. Changing the spiral pattern 



We would like to know how to infer the spiral pattern 
speed from structure in the PBMs. In figure [3] we plot 
the variation of PBMs with a change in the spiral pattern 
speed in the range 0.6 < Sl s /S^o < 1-0; the solar orienta- 
tion with respect to the concave spiral arm is kept fixed 
at <fio — 20°. This range places the Sun from just inside 
the 4:1 ILR to the CR. In contrast to figure [5] here all 
panels have the axisymmetric background subtracted to 
reveal the symmetry of the spiral residuals. We are now 
looking for features that are strong enough to be detected 
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in an actual pencil-beam survey. The z m i n , z ma x values 
indicated above each panel give the maximum error intro- 
duced by spiral structure in the otherwise axisymmetric 
background disk. These values for the line-of-sight ve- 
locity Vd, and its standard deviation a Vd , are ~ 35 km/s 
for vo = 220 km/s (left and middle columns in figure 
[3]), which is well above the resolution of upcoming radial 
velocity surveys (< 3 km/s). Strong features showing 
marked variation with the change in pattern speed are 
apparent in all three observables: 

(1) In the case of Vd high positive and negative velocity 
groups resulting from the effect of the 2:1 ILR are found 
at (l,d/r ) ps (-30°, 0.7) and (l,d/r ) ps (25°, 1.3), re- 
spectively, for Q s — O.6r2o (top left panel). With the 
increase of pattern speed, these clumps spiral in a clock- 
wise direction toward the Galactic center. 

(2) The standard deviation of the line-of-sight velocities, 
a Vd , which can also be described as the " heating" (or the 
random motions) of stars, peaks at a particular ring-like 
shape around the Galactic center for all pattern speeds 
(middle column). These rings are associated with the 2:1 
ILR induced by the spiral density wave. It is clear that 
the radii of these rings are changing with the change of 
the pattern speed and thus the location of the 2:1 ILR. 
Beyond the CR (f2 s = VLq) this resonance falls inside the 
inner three kpc or inside the Galactic bulge. Thus, using 
the radius of this hot ring to infer the location of the 
2:1 ILR (and thus the pattern speed) is only valid if f2 s 
values are in the range considered in figure [3] 

(3) Lastly, the number density PBMs in figure [T] (right 
column) are also indicative of the changing pattern 
speed. The disruption of the spiral arms near the 2:1 
LRs and variation in spiral strength due to the encounter 
the second order resonances and CR, creates a large 
contrast in these axisymmetric background subtracted 
PBMs. Many of these features can be used in addition 
to the information extracted from the velocities. 

All of the features described above can be used to iden- 
tify the location of the 2:1 ILR and thus the patter speed. 
As we mentioned, however, if the solar circle is placed at 
or beyond the CR, the 2:1 ILR falls inside the Galactic 
bulge. Consequently, the features created by it disap- 
pear. Fortunately, just as this happens the 2:1 OLR 
enters the Galactic disk (our disks extend to a radius of 
2ro) and similarly to the 2:1 ILR case, resonant features 
are created, this time in the outer parts of the disk. 

4.2. Changing the Sun's orientation with respect to a 
spiral arm 

How can we infer the Sun's azimuth with respect to 
the galactocentric line passing through the intersection 
of the Solar circle and a concave spiral arm? To find out 
we plot PBMs of simulation runs with the same pattern 
speed and different solar phase angle. Figure [4] shows 
such plots for a fixed Vl s = 0.7^o and a solar phase angle 
changing from top to bottom in the range 0o = [0° , 40°] . 

Similarly to figure [3] we now look for strong features in 
the three observables that can be used to estimate 4>q: 
(1) Inspection of the left column of figure |4] reveals a 
negative velocity stream which changes position with a 
change of phase angle. For <fio — (top left) this fea- 
ture is centered on I ps 0° at a heliocentric distance of 
d/ro « 0.5. As the angle is increased this steam moves 
to larger longitudes roughly preserving its distance from 



the Sun. Note that the high positive and negative fea- 
tures discussed in the context of Vd in figure 02 do not 
vary as the angle is changed since the pattern speed is 
kept fixed. 

(2) The line-of-sight velocity dispersion (middle column 
of figure [3]) does not seem to be particularly useful for 
constraining the solar phase angle. 

(3) Finally, the structure in the number density PBMs 
shows prominent variation with the change in solar angle. 
For example pencil-beam observations at ±45° would be 
drastically different depending on the phase angle. 

In an actual survey we would first try to infer the posi- 
tion of the inner or outer LR as discussed in Section [4~T1 
and thus find the pattern speed. 

5. FOUR- ARMED SPIRAL STRUCTURE 

So far we have only discussed simulations involving a 
two-armed spiral density wave perturbation. In this sec- 
tion we show the effect of a four-armed structure, make 
comparison with the two-armed case and suggest a way 
to distinguish between the two. 

Figure[5]shows PBMs of a galactic disk, similarly to fig- 
ure[31 but perturbed by a four-armed spiral structure. As 
in figure [3] pattern speed changed from top to bottom in 
the range fl s — [0.6, 1.0]f2o and the solar orientation with 
respect to a concave arm is kept fixed at tfio = 20°. In 
this case the first order resonances are the 4:1 ILR/OLR 
which occur at fl s — 0.65r2o, 1.35fio, respectively. 

As expected, more structure is apparent in all PBMs 
in the case of the four-armed structure. Note that the 
hot rings in Acr„ present in the case of the two-armed 
structure (figure [3]) are also apparent in the four-armed 
case although not as pronounced. The reason for this 
is the fact that the 2:1 ILR which causes these is a sec- 
ond order when m = 4. These hot rings can be used in 
both the two- and four-armed cases to estimate £l s but 
are expected to be much stronger for m = 2. Inspecting 
Avd and AE (left and right columns) in figure [3] and \5\ 
it is clear that pencil-beam observations along Galactic 
longitudes I = 0° or I = —90°, for example, can unam- 
biguously distinguish between m = 2 and m = 4 spi- 
ral structure, as the oscillation frequency doubles when 
m = 4. 

6. CONCLUSION 

Upcoming Galactic disk surveys will reveal the age, 
composition and phase space distribution of stars within 
the various Galactic components. These stellar excava- 
tions will provide essential clues for understanding the 
structure, formation and evolution of our Galaxy. To fa- 
cilitate the interpretation of the huge amounts of data 
resulting from these surveys, Galactic disk models, such 
as the one presented here, are needed to interpret the 
observations. 

We have investigated how the Milky Was spiral struc- 
ture parameters, such as pattern speed and solar phase 
angle, can be estimated in a deep all-sky survey. We 
performed a series of test-particle simulations of a warm 
galactic disk approximating the disk kinematics of the 
Milky Way. We considered both two- and four-armed 
spiral structure and suggested a way to distinguish be- 
tween the two using velocity and number density maps. 

We found that the axisymmetric potential needs to be 
known to ~ 10%, line-of-sight velocities to ~ 20 km/s, 
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Fig. 4. — Pencil-beam maps showing variation of structure in the observables with solar orientation at a fixed pattern speed Q s = 0.7f2o- 
As in figure [3] all panels have the axisymmetric background subtracted. The solar phase angle is varied in the range <f>o = [0°,40°]. The 
negative velocity stream found in the Av^ PBMs (left column) at (i,d/rrj) ~ (0°,0.5) moves toward positive longitude with the change of 
phase. This is a strong feature (~ —20 km/s) and can be used to estimate <po. Even more indicative of changes in the solar phase angle 
are the number density PBMs (right column). For example pencil-beam observations at ±45° would be drastically different depending on 
the phase angle. 



and distance uncertainties need to be less than ~ 30%. 
The mean line-of-sight velocity and the velocity disper- 
sion are affected by up to ~ 35 km/s which is well within 
the detectable limit for forthcoming radial velocity sur- 
veys. Pattern speed can be constrained by a hot ring at 
the 2:1 ILR in both two- and four-armed spiral structure. 
To distinguish between the two, however, we also need 
information related to the velocities and stellar number 
density. If the pattern speed is such that the 2:1 ILR is 



hidden inside the Galactic bulge the 2:1 OLR would be 
present in the outer Galaxy and thus can equivalently 
be used to estimate the pattern speed. Once the pattern 
speed is known the solar angle can be estimated from the 
number density variation with heliocentric distance; </>o 
is also reflected in the Vd PBMs. 

Future work needs to address the issue of how to ob- 
tain the axisymmetric background potential needed to 
subtract from the observational data as discussed at the 
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Fig. 5. — Same as figure [3] but in the case of a four-armed spiral structure. Pattern speed changed from top to bottom in the range 
Q s = [0.6, 1.0]f2o and the solar orientation with respect to a concave arm is kept fixed at <j>o = 20°. Inspecting Ai^ and AS (left and right 
columns) it is clear that a pencil beam radial velocity observation along Galactic longitude of £ = —90° can unambiguously distinguish 
between m = 2 and m = 4 structure. 



end of Section |U Also, it is important to know what 
type of tracer stars are needed that would allow the es- 
timation of photometric parallaxes with errors less than 
~ 30%, and the distribution of those stars. 

While here we only considered steady state spiral stric- 
ture, other theories of spiral structure, such as transient 
and swing-amplified spirals, need to be investigated as 
well. It has also been suggested that the Galaxy har- 
bors two sets of s piral structur e moving at the same 
(iLepine et alJl200lD or diffe rent (|Naoz and Shavivll2007l ; 
iMinchev and Q uillen 2006) pattern speeds. We expect 



in all those cases it will be again resonant features to 
relate to the pattern speed and solar angle. In the case 
of non-steady state spirals, however, the structure in the 
PBMs will vary with integration time and interpretation 
will become more complicated. We refer these cases to a 
future study. 

It is also known that the Milky Way is a barred galaxy. 
The simulations performed here do not include the influ- 
ence of the bar. This is not necessarily a shortcoming 
since most of the features in the PBMs we use to in- 
fer spiral structure parameters are caused by resonances 
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and, unless a resonance overlap with the central bar ex- 
ists in the same location, those would not be different 
when a bar is included in the simulations. Future work 
should also look at this problem. 
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